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ABSTRACT 

We suggest a new component of the Milky Way galaxy that can account for both 
the optical depth and the event durations obtained by the MACHO microlensing survey 
toward the Large Magellanic Cloud. This component is consistent with recent evidence 
for a significant population of faint white dwarf stars, detected in a proper motion study 
of the Hubble Deep Field, which cannot be accounted for by stars in the disk or spheroid. 
This new component consists of (mostly) old white dwarf stars distributed in a highly 
extended (very thick) disk configuration. It extends beyond the traditional thin and 
thick disks, but well within the dark, roughly spherical CDM halo. The total mass in 
this component is ~ 7 — 9 x 10 10 M Q . We argue that such a component is reasonable, 
natural, consistent with a variety of observations, and many of the problems associated 
with a significant halo population of white dwarfs are ameliorated. 

Subject headings: dark matter — MACHOs — white dwarfs 
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1. Introduction 

From a comparison of the amount of matter as- 
sociated with the lumino us components of galaxies 
( Faber fc Gallagher 1979) and constrain ts from big 



bang nucleosynthesis jBurlcs et. al 1999] ), it is clear 



that most of the baryons in the Universe today are 
dark. A large fraction of this baryonic dark matter 



may be in the form of h ot, diffuse gas (Gates, Gyuk 



Ho lder fc Turner 199 8r Fukugita, Hogan fc Peebles 



1998), but dark compact objects assumed to reside in 



the halos of galaxies (MACHOs) are also candidates 
for baryonic dark matter. 

The past few years have yielded much exciting 
new data from observational teams searching for evi- 
dence of microlensing along several lines of sight. The 
results from these surveys have raised many ques- 
tions. The microlensing optical depth toward the 
Large Magellanic Cloud (LMC) obtained by t he MA- 
CHO collabor ation, t LM c = (2.1±q\\) x 10 ~ 7 flAlcock 



al. 1997t), is consistent with a significant frac- 

40%) ( |Gates, Gyuk 
The du- 



tion of the galactic halo (20% 
k Turner 1996| ) in the form of MACHOs 
ration of these events also indicates that, under the 
assumption of a spherical isothermal halo distribu- 
tion for the MACHOs, the average MACHO mass is 
~ 0.5M Q , with large statistical uncertainties (Alcock 



2t al. 1997bD - 

Such masses suggest several candidates for the 
lenses including faint halo stars, white dwarfs, and 
black holes. Direct searches for faint halo stars have 
placed severe limits on their contributi on to halo, re- 
quiring it to be less than about 3% (Flynn, Gould 



k Bahcall 1996 



|Graff fc Frcese 1996|) . Primordial 
black hole candidates require a fine tuning of the ini- 
tial density perturbations, and details of QCD phase 
transition black hole formation remain to be worked 



out| in order to assess their viability and mass function 



(Jeciamzik &; Niemeycr 1999). Thus in the standard 



halp model interpretation, wnite dwarls appear to be 
the remaining strong candidate tor tnc lenses. 

However, other possible interpretations of the MA- 
CHO results have been proposed in order to avoid 
the difficulties associated with a large halo popula- 
tion of white dwarfs. The mass estimate for the 
lenses depends upon the assumed phase space dis- 
tribution of MACHOs. There are large uncertainties 
in the halo model parameters, including the distribu- 
tion and velocity structure of the dark matter, and 
attempts have been made to exploit these uncertain- 



ties in order to obtain mass estimates from the current 
data that are consistent with lenses in the substellar 
regime (e.g. brown dwarfs). Previous work by the 
authors and others have examined a wide range of 
halo models, including flattened halos (Gyuk & Gates 
1998| ) , halos with a bulk rotational component to the 



velocity structure ( pyuk fc Gates 1998 ) and halos 
with anisotropic velocity dispersions (Gyuk, Evans & 



Gates 1998). These analyses have shown that for any 



reasonable (smoothly varying) phase space distribu- 
tion of the lenses, the implied lens mass is still much 
larger than the hydrogen burning limit, and thus one 
cannot appeal to modeling uncertainties in order to 
invoke brown dwarfs as candidates for the MACHOs. 
More recent work has found that varying the model 
parameters cannot produce a mass estimate from the 
current data greater than about 0.8A/q, implying 
that neutron stars are also not likely lens candidates 
( Gates, Olinto and Venkatesan 1999 ). 

Other work has explored the possibility that the 
lenses are not in the halo of the Milky Way. LMC 
self-lensing has been suggested by Sahu (1994), but 
recent work by Gyuk, Dalai & Griest (1999) has ar- 
gued that this is unlikely. Zaritsky & Lin (1997) 
and Zhao (1998) have suggested an intervening pop- 
ulation of stars toward the LMC (tidal debris or a 
dwarf galaxy) could be responsible for the microlens- 
ing events. Again this suggestion has been subject to 
much debate, and a recent paper by Gould (1999) ar- 
gues strongly against such a scenario. Galactic mod- 
els in which dark extensions of known populations, 
such as a heavy spheroid or thick disk, could be the 
source of the lenses were explored by Gates, Gyuk, 
Holder & Turner (1998). We will comment further on 
these models in section 3 of this paper. 

However, recent results from a proper motion study 



of the Hubble Deep Field (HDF) (tbata, Richer 



Gilliland & Scott 1999), and from a comparison of 



the n orth and south HDF images ( Mendez & Minnit 



199S ) have added a new piece to the puzzle. These 



studies provide further evidence that there may be a 
previously undetected population of old white dwarfs 
in the galaxy. This strengthens the interpretation of 
the MACHO lenses as white dwarfs, and thus makes 
the above alternatives less appealing. 

Ibata et al.(1999) compared the original HDF with 
a second image of the same field taken approximately 
2 years later, searching for proper motions of faint 
objects. They found 5 faint, blue objects which had 
a significant (> 3<r) shift in the centroid position 
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over this two year period consistent with the detec- 
tion of proper motions of around 20-30 mas/yr. A 
third epoch of observation is planned for approxi- 
mately 2-3 years after the second. Obviously the 
detection of proper motions eliminates the possibil- 
ity that these are extragalactic sources, and indicates 
that they must be relatively close by. 



Using new models for white dwarf cooling (Hansen 
1 9£^9|) , the sources detected by Ibata et al. are con- 
sistent with old (> 12Gyr), 0.5M© white dwarfs. 
Hansen's model predicts that very old white dwarfs 
will be blue and somewhat brighter than earlier mod- 
els in which white dwarfs continue to redden as they 
cool. Previous limits on white dwarfs in the halo were 
based on the older white dwarf cooling models. Ibata 
et al. also argue that these moving sources cannot be 
part of the known disk, thick disk or spheroid pop- 
ulations and further, that the number of sources de- 
tected is consistent with the number expected for an 
all white dwarf halo, although this claim is strongly 
dependent on many factors, especially the assumed 
IMF of the white dwarf progenitor population (Richer 
19S(9|). 

Mendez & Minniti (1999) compared faint blue 
point sources in HDF-North and HDF-South. They 
find 5 such objects in HDF-North and 10 in HDF- 
South. The core of their argument is that this dis- 
tribution is inconsistent with distant extra-galactic 
sources, where an equal number would be expected 
for an isotropic Universe. However, a ratio of ~ 2 is 
roughly consistent with that expected for a galactic 
population since HDF-North looks toward the outer 
Galaxy while HDF-South is pointed more towards the 
center of the Galaxy. Mendez & Minniti also state 
that these sources represent ~ 1/3 — 1/2 of the dark 
matter in the Galaxy. 

While these new data are still somewhat prelim- 
inary, they do raise the intriguing possibility that 
there is a previously undetected population of white 
dwarf stars that are not part of the disk or (known) 
spheroid. These new results, along with the MACHO 
data and the inability of modeling to significantly 
change the mass estimates, seems to be relentlessly 
pointing to white dwarfs as the lenses. So, is the 
halo of our galaxy filled with white dwarfs? A stan- 
dard halo interpretation of these data would say yes 
- a significant fraction of the galactic halo must be in 
white dwarfs. However, such a scenario faces serious 
challenges from many directions, especially given the 
claim that the number of white dwarfs detected by 



proper motion studies is large enough to imply that 
approximately half to essentially all of the halo is in 
the form of white dwarf stars. 

2. White Dwarfs in the Halo? 

When considering the possibility that a large frac- 
tion (or all) of the galactic halo might be in the form 
of white dwarfs, it is extremely important to recall 
the evidence for galactic dark matter, including es- 
timates of the total mass of the Milky Way. A re- 
cent analysis of satellite radial and proper motions 
by Wilkinson and Evans (1999) found a total mass of 
the Galaxy Mtot ~2x 1O 12 M , in good agreement 



with other recent estimates (Kochanek 1996, Zarit 
sky 1998] ) . Wilkinson and Evans also find that the 
halo extends to at least lOOkpc, and possibly much 
further to 150 or 200 kpc. Thus the total mass in 
a white dwarf population that comprises a significant 
fraction of the halo would be of order 10 12 M Q , a num- 
ber which already severely strains the baryon budget 
of the Universe. 

Models which propose such a population must also 
account for the mass in the progenitor population of 
stars and in the metal enriched gas produced during 
the formation of the white dwarfs. Combined with the 
above mass estimate for the galactic halo, such con- 
siderations provide serious challenges for these mod- 
els. For example, consider a white dwarf halo which 
is comprised of at least 50% white dwarfs. The to- 
tal mass in white dwarfs today is thus of order 10 12 . 
The efficiency e(m) for producing a white dwarf from 
a progenitor star of mass m is likely to be 0.25 or 
smaller, depending on the progenitor mass (with an 
upper limit of e = 0.5 for progenitor stars of IMq) 
(Adams & Laughlin 1996). Thus, for a white dwarf 



halo mass of M w d, we expect a mass in the progeni- 
tor population Mp sta rs > ^M wc i and a mass in pro- 
cessed, metal rich gas M zgas > 3M w d- A halo of mass 
Mtot — 2 x 10 12 , half of which is in white dwarfs, re- 
quires a progenitor mass of Mp sta rs > 4 x 10 12 . This 
in turn requires an extremely efficient early burst of 
star formation, through which essentially all of the 
baryons in the Universe are processed. 

^From a cosmological point of view, we can con- 
sider the contribution of the white dwarfs and the 
progenitor population to the matter density of the 
Universe. The Milky Way has a mass to light ratio 
M/L ~ 100 or greater (Zaritsky 1998). If we assume 
that this is a typical value for all galaxies, then galax- 
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ies contribute n g > 100/1200/i = O.O8/1- 1 . Compar- 
ing this with n b h 2 = 0.019 ±0.0024 (95%cl, [Buries ct 



al 1 999), we find that a 50% white dwarf halo exceeds 
the baryon budget ( Slaf | 2 cg ° s ~ 2h) even before con- 
sidering the effects of processing most of the baryons 
through an early star phase. A 20% white dwarf halo 
is also difficult to reconcile with the above estimate 



of since the contributions of the progenitor stars 



will exceed Ob . 

Many authors have explored the implications of a 
halo filled with white dwarfs. These analyses, com- 
bined with the estimates of the total mass in the halo, 
make the possibility of a white dwarf halo even less 
tenable. There are several factors to consider in eval- 
uating such models. 

First, the initial mass function (IMF) of the pro- 
genit or stars must be markedly different than the disk 
IM F (Adams fc Laughlin 1996 , Chabrier, Segretain fc 
Me :a 1996 ). Limits on the IMF arise from both low 
and high mass stars. Low mass stars (< IMq) would 
still be burning hydrogen today and should be visi- 
ble. High mass stars (> 8M Q ) would have evolved 
into Type II supernovae, ejecting heavy metals back 
into the interstellar medium []. From limits on red 
dwarfs in the halo and the galactic metallicity, Adams 
& ^aughlin 1996, Chabrier, Segretain & Mera 1996 
find that the IMF must be sharply peaked about a 
progenitor star mass of in ~ 2M . Adams & Laugh- 
lin conclude that even with the above IMF, the white 
dwarf contribution to the halo is limited to less than 
25% (with 50% being an extreme upper limit). 

Next, the metal enriched gas produced when these 
stars become white dwarfs will pollute the remain- 
ing unprocessed gas, leading to high metallicities 
predicted for the Galactic disk and the interstellar 
medium (into which much of this gas must be blown 
out since the total mass in processed gas is much 



larger than the mass of the disk) (Fields, Matthews 



& $chramm 1997| ). Gibson & Mould (1997) have es- 
timated that the expected amount of C, N and O 
produced would be difficult to reconcile with that in 
pop II white dwarfs. 

The white dwarfs in the halo would also produce 
heavy metals via Type la supernovae. Canal, Isern 
and Ruiz-Lapuente (1997) use this to limit the halo 
fraction in white dwarfs to less than 5 — 10% (or a total 



however, Venaktesan, Olinto & Truran 1999 have argued that 
the bounds on high mass stars are significantly relaxed for pro- 
genitor stars with very low (10~ 4 Zq) or zero metallicity. 



mass in white dwarfs of 5 — 10 x 10 10 Mq). In addition, 
deep galaxy counts limit the fraction of the halo in 
white dwarfs, since the brightly burning progenitor 
stars would be visible ( |Charlot fc Silk ). Finally, it is 
worth mentioning that an all white dwarf halo would 
rule out the existence of other dark matter in the 
Universe (for example cold dark matter) ( pates & 



Turner 1994), opening the door to a host of problems 



with large scale structure formation. 

3. A New Component of the Galaxy 

Given the evidence for a previously undetected 
population of white dwarfs and the severe constraints 
on a halo population consistent with this evidence we 
propose a new component of the Galaxy. Suc h a com- 
ponent was first considered by the authors ( Gyuk & 
Gates 1999 ) in the context of attempting to lo wer the 
mass estimates for the MACHO lenses, and in Gates 



Gyuk, Holder fc Turner 1998 in considering dark ex- 



tensions to known components. 

This new component is essentially a very thick 
(scale height > 2 kpc) population of (mostly) old 
white dwarf stars. It is distinct from known galactic 
populations, both in distribution and age. This "ex- 
tended protodisk" extends beyond the thin and thick 
disk populations, but lies well within the halo. While 
the details of the distribution cannot be determined 
without significantly more data, the general features 
of this proposed model can be illustrated with the 
following example: 

Consider an exponential disk with a volume density 
given by 



p(r, z) 



So 
2h z 



exp{{r - r)/r d )sech 2 (z/h z ) (1) 



where r& = 4.0kpc is the scale length and h z = 2.5kpc 
is the scale height. We assume standard values for the 
position and circular velocity of the Sun, r$ = 8.0kpc 
and v c = 220km/s. 

We also assume a velocity structure, which includes 
a rotational component = 170km/s, of the form 



p(r, <p,z) 



1 



1 1 
M ttti r; 



m ^/(27r) 3 a r a 4> a z 
where the velocity ellipsoid varies as 
<j 2 z = 2nGpohl — nCSoh z . 



(2) 
(3) 
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2a 2 z 



For details on varying these model parameters as well 
as different parameterizations of the generic extended 



protodisk model see Gyuk & Gates f 999 and Gyuk & 



Ga;es 1998. We can also consider a spheroid-like dis- 
tribution for this component. Dynamical estimates 
for the mass of the spheroid are considerably larger 
than the luminous mass, although recent studies of 
the mass function of the spheroid indicate that the 
known spheroid population is unlike ly to be able to 
account for the microlensing events ( |Gould, Flynn fc 



Bal jicall 1998|) . Thus a spheroidal distribution would 
again correspond to a previously undetected compo- 
nent. For such a distribution the total mass is con- 
strained in order not to conflict with the inner rota- 
tion curve of the galaxy, which limits LMC optical 
depths r < 1.2 x 1CT 7 . 

The extended protodisk supports approximately 
half of the local rotation speed, with the remainder 
coming from the thin disk and dark (non-MACHO) 
halo (see e.g. Figure 1). The dark halo in these mod- 
els has a large core radius (> 7kpc) and an asymp- 
totic rotation speed of « 180kpc. The total mass in 
the Galaxy out to 50 kpc is ~ 4.6 x 10 n Af Q . For a 
total mass in the white dwarf extended protodisk of 
M wd = 8 x 10 10 M Q , we find: 

• The optical depth toward the LMC generated 
by this component is r ~ 1.5 x 10 -7 ; 

• The lens mass estimates for the current MA- 
CHO event durations is m ~ 0.4M Q , consistent 
with white dwarf masses; 

• We expect to see roughly twice as many white 
dwarfs in the HDF-North compared to HDF- 
South, similar to the halo models. 

Further, simulations of the proper motions of candi- 
dates in the HDF show results broadly consistent with 
the observations of Ibata et al. (1999). Of course as 
previously stated this is strongly dependent on the 
IMF assumed. 

The main feature of this model, however, is that it 
has a much lower total mass in white dwarfs than halo 
models. As outlined above, it is consistent with both 
the MACHO data and the HDF studies for a total 
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Fig. 1. — Rotation curve for an extended proto 
disk Galactic model. The thick solid line is the 
total rotation curve. Other components are given 
as: dashed=extended protodisk, dotted=thin disk, 
dot dashcd= bulge, long dash=halo. Asymptotic ro- 
tation velocity is ~ 200 km/s, the core radius is 9 kpc, 
the bulge mass is 1.3 x 10 10 M Q . Observational data 
points are taken from Figure 1 of Oiling & Merrifield 
(1998). 

mass in white dwarfs of M w d = 8 x 1O 1O M0. This is 
approximately 1/2 of the mass that would be required 
for a halo distribution of MACHOs which would pro- 
duce the same optical depth. Basically, this reduction 
can be understood because most microlensing is due 
to lenses within about 20 kpc of the Sun for either 
configuration. The extended protodisk has less mass 
beyond that distance than a halo. 

In addition, in these models microlensing takes 
place closer to the observer than in the standard halo 
models and thus where the microlensing tube is nar- 
rower. To obtain the same optical depth the density 
locally must therefore be greater. Thus, in models 
which predict the same microlensing optical depth, 
the number of stars that should be detected in the 
HDF is larger for the extended protodisk models than 
for the halo models. That is, white dwarf counts from 
HDF which imply that 50 — 100% of a standard halo 
is in white dwarfs^ correspond to an optical depth 



2 Of course this interpretation is strongly dependent on the as- 



5 



toward the LMC of r - 3 - 5 x 1(T 7 for the halo 
model. The same number of detected white dwarfs 
corresponds to a lower optical depth (and a smaller 
total mass in white dwarfs) in our fat disk model. 

Finally, the smaller mass in white dwarfs today 
implies a smaller total mass in the progenitor popu- 
lation. For our above example the progenitor mass 
Mp s tars ~ 3.5 x 10 11 M Q , a crucial factor of 10 less 
than that for a 50% white dwarf halo, assuming the 
same IMF in both cases. 

There are several predictions of this model that 
can eventually allow it to be distinguished from a 
standard halo white dwarf population. First, the 
LMC optical depth cannot be much greater than 
about 1.5 x 10~ 7 . Thus if the MACHO and EROS 
observations toward the LMC remain greater than 
2.0 x 10 -7 as the statistics improve, this model would 
be ruled out. Second, because the lenses are con- 
centrated closer to the plane of the galaxy, the typ- 
ical lens-observer distance will be smaller (of order 
5 kpc). This in turn implies an increase in the ex- 
pected number of parallax events (Gyuk & Gates 
19£[9j). Finally, the ratio of optical depths toward 
the Small and Large Magellanic Clouds is expected 
to be of order tsm c/ t lmc ~ 0-8, in con trast with 
tsmc/ t lmc ~ 1.5 ( Sackett fc Gould 1993 ) predicted 
for a standard halo. 



The distribution of event durations ( Gyuk fc Gates 
199^) and a detailed comparison of the distribution 
(in direction and magnitude) of the observed proper 
motions will also differ for a halo vs. extended pro- 
todisk white dwarf population, but these seem un- 
likely to be able to differentiate between models with- 
out significantly more data. 

4. Model Implications 

Because the total mass in the white dwarf popula- 
tion today is significantly lower in this model, some 
of the constraints on a halo white dwarf population 
can be evaded, including those which consider the pro- 
genitor population and the ejected metal enriched gas. 
The total mass in this new component represents only 
about 4% of a total halo mass of 2 x 10 12 . Since essen- 
tially all of the current constraints on white dwarf ha- 
los which limit the halo mass fraction in white dwarfs 
do so at only the 10% level, these constraints can be 



satisfied by our model. This includes the Type la 
supernovae constraints which are dependent on the 
mass in white dwarfs today, and cannot be evaded 
by scenarios which involve somehow hiding the metal 
enriched gas produced by the progenitor stars. 

However, there remains much work to be done to 
more carefully consider the implications of this new 
component. First, we still require an IMF which 
differs significantly from the disk IMF. Assuming a 
log- normal distribution, Adams & Laughlin 1996 and 



sumptions made for the IMF and age of the white dwarf pro- 
genitor population. 



Chabrier, Segretain & Mera 1996 used conservative 
constraints to limit the mass fraction of the high and 
low mass end of the progenitor IMF for a halo white 
dwarf population. While the lower total mass in 
our progenitor population will relax the constraints 
(which are based on the metallicity of the Galactic 
disk) somewhat at the high mass end, the mass frac- 
tion of low mass (m < 1M Q ) stars is constrained by 
number counts of faint low mass stars locally. The 
extended protodisk has an increased local density rel- 
ative to a halo distribution, but a lower total mass, 
resulting in a constraint similar to that for a halo. 
Thus we expect to require a fairly sharp low mass 
drop-off in the progenitor IMF. The implications of 
such an IMF, including the lower fraction of primor- 
dial baryons which is processed through this early 
population, need to be examined in greater depth. 

This new component also provides some intrigu- 
ing hints for cosmology. When did this component 
form and how is it related to galaxy formation sce- 
narios? Can this early starburst population help us 
to trace the baryons in the Universe from their pri- 
mordial state to the present, where we find most of 
the baryons in the intracluster medium? 

5. Conclusions 

We have argued that the microlensing data toward 
the LMC, combined with observations of white dwarf 
stars in a proper motion study of the HDF indicate 
the presence of a new component of the galaxy. This 
component can be generally described as an extended 
distribution that extends at least 2 kpc above the 
galactic plane, but resides well within the halo. It 
is consistent with all data and observations of the 
structure and kinematics of the galaxy, and signif- 
icantly alleviates the considerable problems with a 
halo population of white dwarf stars that is consistent 
with microlensing data. Much work remains to care- 
fully consider the implications of such a component, 
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in particular the formation and evolution of the early 
population of progenitor stars (and resulting metal 
enriched gas) that produced this component. How- 
ever, the significantly lower mass in the progenitor 
population as compared to that for a halo population 
of white dwarfs will allow a reasonable fraction of the 
baryonic mass of the Universe to remain in gas that 
has not been processed through these very early stars. 
Moreover, this component may be a more reasonable 
distribution for the remains of an early starburst pop- 
ulation, in which one would expect a more condensed 
distribution than that of the halo. 
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